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Abstract. 

Using the 8fc x 8A; CCD Mosaic Camera on the Im Schmidt telescope 
in Venezuela, we are conducting a large-scale, deep optical, multiepoch, 
photometric (BVRIHa) survey over ^ 120n° in the Orion OB association, 
aimed at identifying the low mass stellar populations with ages ^ 10 Myr. 

We present initial results for a 34n° area spanning Orion lb, la 
and the B Cloud. Using variability as our main selection criterion we 
derive much cleaner samples than with the usual single-epoch photo- 
metric selection, allowing us to attain a much higher efficiency in follow 
up spectroscopy and resulting in an preliminary list of 74 new low-mass 
(~ 0.4 Mq) pre-main sequence stars. 



1. Introduction 

Crucial aspects of theories of star formation can only be tested by studying 
the stellar populations both in and near molecular clouds. While the earliest 
stages of stellar evolution must be probed with infrared and radio techniques, 
many important questions can only be studied with optical surveys of older 
populations with ages ~ 3-20 Myr. 

In the past, studies of OB associations have been used to investigate se- 
quential star formation and triggering on large scales (e.g., Blaauw 1991 and 
references therein). However, because OB stars are formed essentially on the 
main sequence (e.g., Palla & Stabler 1992) and evolve off it in ~ 10 Myr, they 
are not useful for infering star-forming histories on timescales of 1-3 Myr. More- 
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Figure 1. DSS image of tlie entire 120n° survey area spanning the 
Orion clouds and OB associations, with the molecular complex indi- 
cated by the ^^CO contours from Maddalena et al. (1987). The white 
dashed lines outline the associations, as described in Warren &; Hesser 
(1977). The black dashed lines indicate the strip covered in our first 
observations, for which we report initial results here. 

over, it is not possible to study cluster structure and dispersal or disk evolution 
without studying low-mass stars. Studies of individual clusters in the optical 
and IR have been made (cf. Lada 1992), but these are biased toward the highest- 
density regions, and cannot address older and/or widely dispersed populations. 

Recent technological advances have now made it possible to carry out large- 
scale studies, building on the availability of cameras with multiple CCDs on 
telescopes with wide fields of view. 

Figure 1 is shows the Orion A and B clouds and surroundings which we 
propose to survey. The Orion belt stars, 6, e, and are shown for reference. 
The prominent bright emission nebulae are the Orion Nebula (ONC), NGC 2023, 
and NGC 2024 clusters, marking the sites of very recent star formation. Also 
indicated are the OB associations in the region (Blaauw 1964): Ori lb and Ori 
la. Ori Id corresponds to the Trapezium/ONC region; Ori Ic is the region 
surrounding it. Photometric analysis of the O, B and A stars (Warren &; Hesser 
1977, 1978; Brown et al 1994, BGZ) indicate ages of < 1 Myr (Id, see also 
Hillenbrand 1997), 3 Myr (Ic), 7 Myr (lb), and 12 Myr (la). The latest results 
from Hipparcos (de Zeeuw et al 1999) yield mean parallaxes corresponding to 
330 pc (la), 440 pc (lb), and 460 pc (Ic), with an uncertainty of 30%. 
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Figure 2. Left: Diagram of the QuEST Camera finger assembly. 
Each "finger" of 4 CCDs has its own filter. Stars enter the array 
through finger 4 and exit on finger 1. Right: The 1 a dispersion in 
the V magnitudes in a QuEST dataset for 16 scans in Orion. Crosses 
indicate variable stars (99.99% confidence). Most stars 96%) are 
non- variable (dots). 

2. The Photometric Variability Survey 

The large scale, multiband (BVRIHa), multi-epoch, deep photometric survey 
over 120n° in Orion (Figure 1) is being carried out using the QuESTQ camera, 
a 8k X 8k CCD mosaic detector installed on the Im (clear aperture) Schmidt 
telescope at Llano del Hato Observatory, in the Venezuelan Andes (8°47' N, 
3610 m elevation). The 16 2048 x 2048 UV-enhanced, front illuminated, Loral 
CCD chips are set in a 4x4 array (Figure 2, left), covering most of the focal 
plane of the Schmidt telescope, yielding a scale of 1.02" per pixel and a field 
of view of 2.3° x 2.3°. The camera is optimized for drift-scan observing in the 
range —6° < 6 < +6°: the telescope is fixed and the CCDs are read out E-W at 
the sidereal rate as stars drift across the device, crossing each of the four filters 
in succession. This procedure generates a continuous strip (or "scan") of the 
sky, 2.3° wide; conversely, one can survey the sky at a rate of 34. 50° /hr / filter , 
down to Viirn = 19.7 {S/N = 10). 

2.1. Data reduction and analysis 

QuEST has developed its own software since the huge amount of data produced 
is very difficult to handle with packages such as IRAF. The whole process is 
completely automated with minimum interaction from the user. The output 
catalogs contain, among others, J2000.0 positions, instrumental and calibrated 
magnitudes in 4 bands, and their corresponding errors. Positions are good to 
±0.2" down to V-^ 19, within a few square degrees. 



^ The QuEST (Quasar Equatorial Survey Team) collaboration includes Yale University, In- 
diana University, Centro de Investigaciones de Astronomia, and Universidad de Los Andes 
(Venezuela). The main goal of QuEST is to perform a large scale survey of quasars. 
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Figure 3. V vs. V-Ic diagram for different samples, (a) Upper left: 
all stars in Ori lb; (b) upper right: all variables in Ori lb; (c) lower 
left: all stars in a field far from the clouds; (d) lower right: all variables 
in the same field. Solid line: ZAMS. PMS variables in lb are clearly 
separated from the background, but they do not show in the control 
fields. 



We have developed tools for identifying variable stars using differential pho- 
tometry. Using a test and assuming a Gaussian distribution for the errors, 
we consider variable only those objects for which the probability that the ob- 
served distribution is a result of the random errors is very small. In Figure 2, 
right, we show a sample result from the variability analysis in the 2.3" wide strip 
indicated in Figure 1. The dispersion increases for fainter objects, so that most 
(non-variable) objects populate a curved region. We use the test to identify 
potential variable stars with a 99.99% confidence level (crosses in Figure 2). 
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Figure 4. FAST spectra for 2 new MO Orion T Tauri stars (TTS). 
(a) Upper left: Classical T Tauri star (CTTS - W[Hq;]= -32A). Other 
lines arc also seen clearly in emission (H/3, [O I], etc), (b) Upper right: 
the region around 6700A. The Li IA6707 and Ca IA6717 absorption lines 
can be clearly seen. The prominent emission is He IA6678. (c) Lower 
left: Weak line T Tauri star ( WTTS - W[Ha]= -L7A). (d) Lower 
right: region around 6700A, showing strong Li IA6707 aobsorption. 



3. Results 

During Dec. 98 - early Feb. 99, we obtained 16 BVRI scans over the strip indicated 
in Figure 1. We calibrated a 'master' scan using Landolt (1992) standard fields 
and then normalized the photometry in all the other scans to this reference scan. 

The value of variability for picking out pre-main sequence (PMS) candidates 
is shown in Figure 3. The upper left panel is a color-magnitude diagram with 
all the stars in a ~ 10n° area within Ori lb (the ZAMS at 440 pc is shown for 
reference); the upper right panel shows the variables in the same field, picked out 
using our selection criteria. Populations above and below the ZAMS are clearly 
separated using our variability criteria. Indeed, the only 5 known TTS in this 
region (Hcrbig & Bell 1988) were all recovered as variables. We also compared 
our data with the Kiso Ha survey (c.f. Wiramihardja et al. 1993) and found 
that the vast majority of Ha stars above the ZAMS are detected as variables 
(-^ 70%), but essentially none of the ones below are (they could be a mixture of 
field dMe stars and some false detections; c.f. Briceno et al. 1999). 

To emphasize the point even further, in the lower panels of Figure 3 we 
show color-magnitude diagrams for field stars {a = 4h — 5h, 6 = —1°), showing 
clearly that the tail of the distribution of background stars extends far above the 
main-sequence, while variables show essentially no population of PMS objects. 
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We have initiated followup spectroscopy of the brighter (V ^ 16.5) variable 
PMS candidates (Fig.3b), using the FAST spectrograph (Fabricant et al. 1998) 
on the 1.5m telescope at SAO, with a spectral resolution of 6A in the range 4000 
- 7000 A. We confirm low mass PMS stars based on the presence of emission in 
Ha and other lines, and of Li IA6707A strongly in absorption, which is a reliable 
indicator of youth in K4-K5 and later spectral type stars (c.f. Bricerio et al. 
1997). Even at this low resolution, Li I can be seen in late type stars with 
high SNR spectra (Figure 4). In this way, we have obtained spectra for 157 
candidates and confirmed 74 of them as new TTS. We are now placing these 
objects in the HR diagram to derive their masses and ages. This high (~ 50%) 
efficiency is the result of the clean selection provided by the variability criterion. 

The new TTS have spectral types K7 - M2, corresponding to masses roughly 
0.5 - 0.3 Mq at ~ 1 - 3 Myr. Though preliminary, this list of new TTS already 
suggests that the fraction of CTTS in la is much lower than in lb, which would 
be expected if la is indeed older than lb. We are analyzing in detail the light 
curves of these new stars, and spectroscopy of further candidates is under way. 

This is the first optical survey to approach in spatial extent studies of ex- 
tended star-forming regions (near the galactic plane, not reached by SDSS) done 
by other surveys like the RASS and 2MASS, but going much deeper than the 
RASS and having simultaneous photometry over several optical bandpasses for 
many epochs, providing a unique variability database that one-time surveys like 
2MASS cannot not offer. 
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